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Abstract 

In the paper we discuss possible applications of the so-called stroboscopic tomography (stro¬ 
boscopic observability) to selected decoherence models of 2-level quantum systems. The main 
assumption behind our reasoning claims that the time evolution of the analyzed system is given 
by a master equation of the form p = L p and the macroscopic information about the system 
is provided by the mean values mj(tj) = Tr(Qip(tj )) of certain observables {Qi}I=i measured 
at different time instants {fjg =1 . The goal of the stroboscopic tomography is to establish the 
optimal criteria for observability of a quantum system, i.e. minimal value of r and p as well as 
the properties of the observables {Qi}l =1 . 


1 Introduction 

According to one of the most fundamental assumptions of quantum theory, the density matrix car¬ 
ries the achievable information about the quantum state of a physical system. In recent years the 
determination of the trajectory of the state based on the results of measurements has gained new 
relevance because the ability to create, control and manipulate quantum states has found applica¬ 
tions in other areas of science, such as: quantum information theory, quantum communication and 
computing. 

The identification of an unknown state by appropriate measurements is possible only if we have 
a set of identical copies of this state, because each state can be measured only once due to the fact 
that every measurement, in general, changes the state. Moreover, in order to create a successful 
model of quantum tomography one needs to find a collection of observables, such that their mean 
values provide the complete information about the state. In the standard approach to quantum 
tomography of 2-level systems one takes as the observables the set of Pauli matrices, denoted by 
{(Ti, < 72 , 03 }, see for example [1]. The reconstruction of the initial density matrix is possible due to 
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the decomposition in the basis {I,<ri, 02 , 03 }, which has the form 



where s t is the expectation value of 0 * in state p( 0). Thus in the standard approach one needs to 
measure three different physical quantities in order to reconstruct the density matrix of a 2 -level 
system. In general for an N-level system one would need to measure N 2 — 1 different observables - 
more about general approach can be found in [2, 3]. This fact implies that the standard approach 
seems rather impracticable as from experimental point of view it is difficult to find as many different 
observables. 

Therefore, in this paper we follow the stroboscopic approach to quantum tomography which 
was proposed in [4] and then developed in [5, 6 ]. In the stroboscopic approach we consider a set of 
observables {Qi}i = i (where r < N 2 — 1) and each of them can be measured at time instants {tj}^ =1 . 
Every measurement provides a result that shall be denoted by mi(tj ) and can be represented as 
rrii(tj ) = Tr(Qip(tj)). Because in this approach the measurements are performed at different time 
instants, it is necessary to assume that the knowledge about the character of evolution is available, 
e.g. the Kossakowski-Lindblad master equation ["] is known or, equivalently, the collection of Kraus 
operators. Knowledge about the evolution makes it possible to determine not only the initial density 
matrix but also the complete trajectory of the state. To make this issue clearer from now on we 
assume the following definition [ 6 ]. 

Definition 1.1. An N-level open quantum system is said to be (Qi, ...,Q r )-reconstructible on an 
interval [0, T] if there exists at least one set of time instants {tj}^ =1 ordered as 0 < t\ < ... < t p < T 
such that the trajectory of the state can be uniquely determined by the correspondence 

[0, T] 3 tj -3 mi(tj ) = Tr(Qip(tj )) ( 2 ) 


for i = 1,..., r and j = 1, ...,p. 

The outcomes that we obtain from the measurements can be presented in a matrix form as 



mi(t 2 ) •• 

• rni(t p ) 

m 2 (ti) 

m 2 (t 2 ) ■■ 

.. m 2 {t p ) 

m r (ti) 

m r (t 2 ) •• 

■ m r {t p ) 


The fundamental question that we formulate is: Can we reconstruct the initial density matrix p(0) 
for a given master equation from the set of measurement results presented in (3) ? 

Other questions that arise in this approach concern: the minimal number of observables for a 
given master equation and their properties as well as the minimal number of time instants and their 
choice. The general conditions for observability have been determined and will be presented here 
as theorems and the proofs can be found in papers [4, 5, 6 ]. 


Theorem 1.1. For a quantum system which evolution is given by the Kossakowski-Lindblad master 
equation of the form 



(4) 


where the operator L is called the generator of evolution, there exists a number (denoted by p) which 
is called the index of cyclicity and is interpreted as the minimal number of observables required to 
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reconstruct the density matrix. The index of cyclicity can be computed from the equality [4] 

r) := max {dimKer( L — AH)}, (5) 

Aeo-(L) 

where cr(L) denotes the spectrum of the generator of evolution (i.e. the set of all eigenvalues of h). 

According to theorem 1 for every generator of evolution there always exists a set of r) observables 
such that the system is (Qi, ...Q^j-reconstructible. Moreover if the system is also (Qi, 
reconstructible, then rf > p. The index of cyclicity seems the most important factor when one 
is considering the usefulness of the stroboscopic approach to quantum tomography. This figure 
indicates how many distinct experimental setups one would have to prepare to reconstruct the 
initial density matrix in an experiment. The index of cyclicity is a natural number from the set 
{1,2,... ,N 2 — 1} (where N = dimR) and the lower the number the more advantageous it is to 
employ the stroboscopic approach instead of the standard tomography. Moreover one can notice 
that the index of cyclicity can be understood as the greatest geometric multiplicity of eigenvalues 
of the generator of evolution. Thus, one can conclude that the index of cyclicity has physical 
interpretation, which is important from experimental point of view, but its value depends on the 
algebraic properties of the generator of evolution. Therefore, the question whether the stroboscopic 
tomography is worth employing or not depends primarily on the character of evolution of the 
quantum system. 

Another problem that we are interested in relates to the necessary condition that the observ¬ 
ables ( Qi,...Q r ) need to fulfill so that the system with dynamics given by (4) will be (Qi, ■■■Q r )- 
reconstructible. First, we introduce a denotation B(fH ) which shall relate to the vector space of all 
linear operators on R. Then by ( A\B) we shall denote the inner product in this space, which is 
defined as 

(A\B) = Tr(A*B). (6) 

Furthermore, one can notice that assuming the dynamics given by (4) the formula for p{t) at 
an arbitrary time instant can be expressed in terms of a semigroup 


m- i 

p{t) = exp(Lt)p(0) = ^a fc (t)L fc p(0), (7) 

k =o 

where p stands for the degree of the minimal polynomial of L. 

This observation enables us to expand the formula for results of measurements (see (2)) in the 
following way 


fi—1 [l— 1 

miitj) = ( Qi\p(tj )) = ^2a k {t j )(Q i \L k p(0)) = ^ a k{tj){(^*) k Qi\p(0)), ( 8 ) 

k =0 k =0 

where L* is the dual operator to L or, in other words, L in the Heisenberg representation. 

It can be proved that the functions afc(i) are mutually linearly independent and can be computed 
from a system of differential equations [ 6 ]. Therefore, the data provided by the experiment allows 
us to calculate the projections ((L*) fc Q,;|/ 9 ( 0 )) for k = 0,1, • ■ ■ , p — 1 and i = 1,2, • ■ ■ , r. It can be 
observed that the initial state p( 0 ) (and consequently the trajectory exp(ht) p(0)) can be uniquely 
determined if and only if the operators (L *) k Qi span the vector space of all self-adjoint operators 
on R. This space shall be denoted by B*{R) and will be referred to as the Hilbert-Schmidt space. 
Now, if the evolution of the system is given by (4) the conclusion can be presented as a formal 
theorem. 
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Theorem 1.2. The quantum system is (Qi, ...Q r )-reconstructible if and only if the operators {Q \,..., Q r } 
fulfill the condition [4, 5] 

r 

Q)K fl (h,Q i ) = B*(H), (9) 

i =0 

where ® denotes the Minkowski sum of subspaces, p is the degree of the minimal polynomial of L 
and K^{L,Qi) denotes Krylov subspace, which is defined as 

K^,Qi) :=Span{Q i ,lAQ i ,{lA) 2 Q i ,...,(h*y- 1 Q i }. (10) 

Remark 1.1. In the theorem 2 we denote by Qo an identity matrix of the appropriate dimension. 

One can notice that for any generator of evolution L we have K^fL,!) = I. 

When discussing the usefulness of the stroboscopic tomography, it can be observed that if one 
takes a hermitian operator Q which belongs to the invariant subspace of the Heisenberg generator L*, 
then K^(L, Q) = Q. Therefore, multiple measurement of the same observable leads to projections 
of /?( 0 ) into distinct operators only if the observable does not belong to the invariant subspace of 
L*. Thus if one considers the implementation of the stroboscopic tomography in an experiment, its 
effectiveness depends on whether one can measure such a quantity that the corresponding hermitian 
operator does not belong to the invariant subspace of the Heisenberg generator. 

The last theorem which will be presented in this section gives the condition for the choice of 
time instants. 

Theorem 1.3. The determination of the initial state of the quantum system with evolution given 
by (4) and which is (Qi, ...Q r )-reconstructible is possible if the time instants {tj} l - =1 satisfy the 
condition [6] 

det [a k (tj)] 0 , ( 11 ) 

where k = 0,1,..., fi — 1. In the above relation ak{tj) denotes the functions that appear in the 
polynomial representation of the semigroup <f>(t) = expfht). 

Having summarized the most important general results concerning the stroboscopic tomogra¬ 
phy, we can proceed to analyzing specific examples. In the main part of this article there are three 
different decoherence models of 2 -level quantum systems, to which the stroboscopic approach has 
been applied. Section 2 is devoted to the problem of dephasing, in which the stroboscopic approach 
allows us to give the concrete formula for the initial density matrix. In section 3 we discuss the 
usefulness of the stroboscopic approach in case of depolarization, which is another model of deco¬ 
herence. Finally, in section 4 we tackle a more general problem, where the stroboscopic approach 
seems to have the greatest advantage. In that section we introduce a parametric-dependent family 
of Kraus operators for which the generator of evolution has no degenerate eigenvalues, i.e. in that 
case there exists one observable the measurement of which performed at three different instants is 
sufficient to reconstruct the initial density matrix. 

2 Quantum tomography model for dephasing 

In this section we are analyzing a decoherence model which is called dephasing. It is a model of 
two-level atoms subject to fluctuating external magnetic or laser fields. In geometric language it 
refers to shrinking of the Bloch ball in x and y directions, and 2 : being left intact. 

The canonical Kraus operators have the following forms [8]: 

K 0 (t) = K^t) = 3, ( 12 ) 
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where n(t) is a function which depends on time and can be expressed as n(t) = e -7 *, where 7 G M+ 
is a dephasing parameter. 

The collection of Kraus operators constitutes a completely positive and trace-preserving map 
(i.e. it is a quantum channel). Therefore, p{t) at any time can be obtained from the formula 


Pit) =^2 K i{t)p{0)K*(t) = 1 + K ^ p( 0) + 0)o~3- 


(13) 


i=0 


Having the specific form of Kraus operators one can obtain the Kossakowski-Lindblad equation 
for evolution of such a system by differentiating the equation (13). It has the following form: 




(14) 


Here, since 03 is a self-adjoint operator, it can also be presented in terms of a double commutator 

(15) 




The explicit form of the generator of evolution can be obtained by using the relation from 
vectorization theory [9] 

vec(ABC) = (C T < 8 > A)vecB , (16) 

where it is assumed that the matrices A,B,C are selected in such a way that the matrix product 
ABC is computable, i.e. the corresponding sizes of A,B,C are si x S 2 ,S 2 x S 3 and S 3 x s 4 , where 
Sl,S2,S3,S 4 e N. 

Taking into account this property one can get the generator of evolution for this system in the 
matrix form 

0 0 0 0 ' 

0 10 0 
0 0 10 ' 

0 0 0 0 


L = —7 


(17) 


Calculating the characteristic polynomial of this operator one obtains 

detIJL — XI) = A 2 (A + 7 ) 2 , 


(18) 


from which one can observe that the index of cyclicity of the system with such a generator is equal 
2 . According to the theorem 1.1 it means that there exist two observables the mean values of which 
enable us to reconstruct the initial density operator and, as a result, the whole trajectory of the 
state. We can instantly notice the benefit of the stroboscopic approach in comparison with the 
standard model - here we measure two different quantities instead of three. 

Furthermore, it can be noticed that operator L fulfills an equality 

L 2 + 7 L = 0, (19) 

which means that deg /i(A,L) = 2, where by p(A, L) one should understand the minimal polynomial 
of L. 

The observables Qi,Q 2 that are needed to perform quantum tomography have to satisfy the 
necessary condition, which is 

2 

QK 2 (L,Q i ) = B*(U) i ( 20 ) 

i=0 
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where K 2 (L, Qi ) denotes a Krylow subspace and can be rewritten as 

K 2 (L,Q i ) = Span{Q i ,L*Q i }. 


( 21 ) 


If we take Q\ = a\ and Q 2 = o 2 + < 73 , it can be observed that h*Q\ = —jQi = —jai and 
L*Q 2 = —ycT 2 . Then as I, Qi,Q 2 and L*Q 2 are linearly independent the condition (20) is fulfilled. 

Therefore, in order to reconstruct the density matrix p( 0) for the system which evolution is 
given by (15) it is enough to measure Q\ once and Q 2 twice for different time instants t\ and t 2 . 

To obtain a specific formula for the density matrix we have to analyze the polynomial rep¬ 
resentation of the completely positive map 4>(t) = exp(Lt). It has already been mentioned that 
deg /i(A, L) = 2 , thus the polynomial representation takes form 

exp(ht) = ao(t)I + ai(t)L. (22) 


Coefficients «j(t) that appear in this equation can be easily computed because they satisfy a set of 
differential equations [ 6 ], which has been mentioned in the introductory section. One can easily get 


ao = 1 and a\{t) 


-( 1 -e-*). 
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(23) 


Bearing in mind the general representation of a result of measurement (see equation ( 8 )) we 
shall write the formula for m 2 (t\) 


m 2 (ti) = ( Q 2 \p(h )) = (Q 2 \exp(hti) p(0)) = ^ a fc (ti)((L*) A: Q 2 |p( 0 )). (24) 

k =0 


As the observable Q 2 is going to be measured twice, we obtain two similar equations for the two 
results. They can be combined into a matrix equation 


m 2 {ti) 


_m 2 (t 2 ) 



1 


1 


^( 1 

1 i(l 



(QMo)) ' 

(l*q 2 |p(o))J ’ 


(25) 


where m 2 (t\) and m 2 (t 2 ) are results obtained from measurement, i.e. mean values of the observable 
Q 2 in two different time instants: m 2 {ti) = Tr(Q 2 p(t l )). It is evident that if t\ / t 2 , 


det 


1 i(l-e-^) 
1 i(l-e-7fe) 


/ 0 . 


(26) 


Therefore, from equation (25) we can calculate projections of p(0) into the operators Q 2 and h*Q 2 . 
We obtain the following results 


(QMo)) 


m 2 (ti){l-e ' yt2 ) - m 2 (t 2 )(l - e ltl ) 
e -ih — e ~T2 


(27) 


<L*Q 2 |p(0)) = 


7(m, 2 (t 2 ) - m 2 (t 1 )) 


(28) 


— e ~T2 

For the observable Q 1 it is sufficient to write one equation and transform it in the appropriate 


way 


= (Qi|p(0)) + -(1 - e 7tl )(L*Qi|p(0)) = 

7 

= <Qi|p(0)> + ^(1 - e-^)(- 7 )(Qi|p(0)) = (29) 

= ( 1 + ( e -^-l )) (Q 1 |p(0 )) = 

= e-^(Qi|p(0)), 
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from which we get the projection of p{ 0) into operator Q\ 

(Qi\p{0)) = m 1 (ti)e ytl . (30) 

In order to obtain an explicit formula for the density matrix we will use the theorem that any 
two-dimensional density matrix can be expanded using the identity operator I and the traceless 
Pauli matrices { 07 , <r 2 , < 73 }. The decomposition takes the from 

p{ 0 ) = ^ ’ ( 31 ) 

where s* = Tr(aip( 0 )) = (cjj|p( 0 )). 

We can notice that 


(Q 2 |p( 0 )) = (02 + < 7 3 |p( 0 )) = <o- 2 |/o(0)> + (<r 3 |p(0)), (32) 

(L*Q 2 |p(0)) = —7(cr 2 |p(0)). (33) 

Taking these equations into account we get the projections we need to reconstruct the density matrix 

m 2 (t 1 ) - m 2 (t 2 ) 


(o- 2 |p(0)) = 


=-7*l — g~7*2 


and 


<cr 3 |p(0)> = 


m 2 (t 2 )e 


-7*1 _ 


m 2 (ti)e 


-7*2 


(34) 


(35) 


e -7*i _ e -7*2 

Having found all the projections of p(0) one can write an explicit formula for the density matrix 


1 


p{ 0 ) = ^ ( I + mi(fi)e 7tl cri + 


m 2 (t\) — m 2 (t 2 ) m 2 (t 2 )e 


-7*1 _ 


7*1 — p ~ 7*2 


o 2 + 


m 2 (ti)e 


-7*2 


=—7*1 — p — 7*2 




(36) 


which is the final result of this analysis. 

In this section we were capable of creating a successful model of quantum tomography by using 
the stroboscopic approach. The advantage over the standard model is the fact that here one needs 
to measure only two different observables instead of three. The number of pairs of observables that 
fulfill the necessary conditions of observability for the system in question is infinite. Nevertheless, 
it has been shown for the selected two observables that it is possible to obtain the explicit formula 
for the initial density operator. 


3 An attempt to apply the stroboscopic tomography to depolar¬ 
ization 

Depolarization is another model of decoherence that is being analyzed in this paper. Geometrically 
speaking, this model refers to squeezing the Bloch ball uniformly in the three directions. The Kraus 
operators for depolarization have the forms [8] 

Ko(t) = \J 1 + { m) = , (37) 

l J i=l, 2 ,3 

where n(t) depends on time and can be expressed as n(t) = e _7< , where 7 € M_|_ is called depolarizing 
parameter. 


7 











A. CzerwihsM ONE-PARAMETRIC NON-DEGENERATE FAMILY OF KRAUS OPERATORS 


Having the collection of Kraus operators for this case one can calculate the derivative of p 
analogously as in section 2. It leads to the evolution equation in the Kossakowski-Lindblad form. 
One can obtain 

^ = ai P ai - 3 ^) • ( 38 ) 

As in this case the operators {<7j}i=i, 2 ,3 that govern the evolution are hermitian, the master equation 
can be represented as a sum of double commutators 


dp 

dt 


“g 

i— 1 


(39) 


Applying the relation from vectorization theory (16), one can get the explicit form of the gen 
erator of evolution 

7 

4 

which can also be presented in the matrix form 


L = j (ct 3 ® 07 + cnf (g> <j 2 + cr 3 <g> cr 3 — 3I 4 ) 


(40) 


L = — — 


10 0-1 
0 2 0 0 

0 0 2 0 

-10 0 1 


Now one can find the eigenvalues of the operator L 

Ai = —7 and A 2 = 0, 

and their corresponding multiplicities which are 

ni = 3 and n 2 = 1. 


(41) 


(42) 

(43) 


One can quickly check that there are three linearly independent eigenvectors of L that corre¬ 
spond to the eigenvalue A 3 . It means that the index of cyclicity for the generator of evolution given 
by (41) is equal 3, which implies that we need 3 different observables to perform quantum tomogra¬ 
phy on the system. Therefore, in case of depolarization stroboscopic approach to tomography has 
no advantage over the standard tomography model for a 2-level system. 


4 One-parametric non-degenerate family of Kraus operators 

In this section, before we introduce the main result, we shall assume two definitions. One concerns 
a collection of Kraus operators that can be associated with any completely positive map. The other 
one relates to a family of Kraus operators which should be understood as parametric-dependent 
and, therefore, more general case. To illustrate the difference between these two terms we shall first 
revise the theorem on completely positive maps [10]. 

Theorem 4.1. A linear map 4>(f) : B(fH) —> B(TL) is completely positive if and only if it is of the 
form 

C 

m(X)=^K i (t)XK*(t), (44) 

i=l 


where Ki(t) G B(TL). 
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Now we can make a distinction between a collection of Kraus operators and a family of Kraus 
operators. 

Definition 4.1. A set of operators {Kflt )}^ =1 that appear in the above decomposition of a completely 
positive map is called a collection of Kraus operators. 

Definition 4.2. If we have a completely positive map <£(£,0), where t stands for time and a £< 
0"min^CL ma x > is a real parameter that influences the structure of the generator of evolution L asso¬ 
ciated with the map, and we denote the sum representation o/<b(t,a) as 

C 

4>(f, a)(X) = 2 Kflt, a)XK*(t, a), (45) 

i=1 

then we shall call the set of operators {K l (t, a)}^ =1 a one-parametric family of Kraus operators on 
the interval <C Umm, a max P • 


Remark 4.1. According to the definition of the one-parametric family of Kraus operators, the 
parameter has to influence the structure of the generator L and, consequently, the positive constant 
7 from the two previous sections does not create a family of Kraus operators as the generator L is 
proportional do 7 (see (17) and (41)j. 

Remark 4.2. In a similar way we could introduce the definition of k-parametric family of Kraus 
operators for k = 2, 3, - - - 


One can notice that according to these definitions a family of Kraus operators comprises an 
infinite number of collections of Kraus operators. 

Now we can proceed to the main part of this section. As the third model of decoherence we are 
analyzing the case when the evolution of the open quantum system is given by a one-parametric 
family of Kraus operators 

/1 + 2 n(t) . . / all — n(t)) 

K 0 (t,a) = y---1, K\(t,a) = y--- -04 

and K 2 (t,a) = 

where n(t) depends on time according to n(t) = e _7t , where 7 G M + is a decoherence parameter, 
and a is a real parameter that influences the structure of the generator of evolution. 

Let us observe that it is a family of Kraus operators for a£l, Nevertheless, some constraints 
need to be found concerning a because the family of Kraus operators from (46) have to constitute 
a completely positive map which is strictly trace-preserving as the evolution should retain all the 
properties of the density operator. One can notice that the sufficient and necessary conditions for 
this are 

a > 0 and 2 — a > 0. (47) 

It is easy to check that under these two condition the following equality holds 


(2 - o)(l - /e(i)) 


- 02 , 


2 


i =0 


(48) 


thus any collection of Kraus operators taken from the family (46) with a £< 0; 2 > constitutes a 
quantum channel and can be treated as a model of decoherence. 

In this section we propose the following theorem concerning the introduced family of Kraus 
operators. 
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Theorem 4.2. If the evolution of an open quantum, system is given by a collection of Kraus oper¬ 
ators from the family defined in (46) for any a G (0; 2) \{1}, then the eigenvalues of the generator 
of evolution L are non-degenerate. 

Proof. Kraus operators allow us to find the equation for evolution of the system, which takes the 
following form 

dp 2 a 2 — a 

= -g7P+ H-— 7 a 2 pcr 2 . 

Using again the idea of vectorization (16) the explicit form of the generator L can be written as 

L = ^ (— 2 I 4 + aai <S> o\ + (2 — a)a\ ® 02 ) . (50) 

Consequently its spectrum can be computed 

o-(L) = {-^ 7 , 0 , ^( 0 - 2 ) 7 ,-^ 07 }. (51) 

Now it can be observed that for any a G (0; 2) \{1} the spectrum of the generator L consists of 
four different eigenvalues, which completes the proof. □ 

The theorem presented in this section shows that there exists a one-parametric family of Kraus 
operators for which eigenvalues of the generator of evolution L are non-degenerate, which means 
that the index of cyclicity is equal 1. Henceforth a family that possesses this property shall be 
called a one-parametric non-degenerate family of Kraus operators. 

The result means that for any evolution given by (46) with a G (0; 2) \{1} there exists one 
observable the measurement of which performed at three different instants allows us to reconstruct 
the initial density matrix and, as a result, the trajectory of the state. From experimental point of 
view the stroboscopic approach seems to have a considerable advantage over the standard tomog¬ 
raphy as in this case an experimentalist needs to prepare only one kind of measurement and repeat 
it three times instead of performing three different measurements. 

5 Summary 

This paper gives a brief insight into the possible applications of the stroboscopic tomography to 
2-level decoherence models. It has been shown that the usefulness of this approach differs depending 
on the kind of evolution. The most promising result is described in section 4. In that part it has been 
proved that there exists a one-parametric non-degenerate family of Kraus operators, which means 
that for an infinite number of generators of evolution the index of cyclicity is equal 1 , i.e. for any 
a G (0; 2) \{1} there exists one observable the measurement of which performed at three different 
time instants provides sufficient data to determine the trajectory of the state with dynamics given 
by (46). Further research into the problem of non-degenerate families of Kraus operators is planned 
for the foreseeable future. 
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